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1 Introduction 

> ' 

Q\ ; A good knowledge of total crosss-sections of high energy photon-proton and photon- 

photon interactions, in an energy domain where no data are available, is important 
^j- ■ from the point of view of understanding high energy cosmic ray data and planning of 

, \ ', the TeV energy e + e~ colliders respectively. This requires developing a model which 

can explain the current energy dependence of these cross-sections observed in the 
laboratory experiments and then using it to predict the cross-sections in the required 
higher energy regime. Apart from this very prosaic reason for studying the subject, 
high energy behaviour of total hadronic cross-sections is an issue of great theoretical 
importance. Very general arguments based on unitarity, analyticity and factorisation 
in fact imply a bound on the high energy behaviour of total hadronic cross-sections [1] . 
This bound predicts, independent of the details of the strong interaction dynamics, 
that asymptotically a to t < C(logs) 2 . All the experimentally measured hadronic cross- 
sections seem to rise with energy [2], although it is not clear whether the rate is the 
same for all the hadronic processes; nor is it clear whether the asymptotic behaviour 
is already reached at the current energies. In view of the important clues to the strong 
interaction dynamics that this energy dependence holds and the equally strong need of 
its precise knowledge in the high energy regime for the planning of future experiments, 
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or the understanding of the high energy cosmic ray data, it is not surprising that this 
has been a subject of intense theoretical investigations [3j 0]. 

Since QCD has now been established as the theory of strong interactions, it is of 
course of import to seek an understanding of this interesting bound in a QCD based 
picture and try to see clearly as to which piece of the QCD dynamics is it most closely 
related to. A description of one effort j5] in this direction, viz. a QCD based model to 
describe the energy dependence of the total hadronic cross-section and its extension [6] 
to the case of the high energy behaviour of the photon induced processes is the subject 
of this note. We will first summarise very briefly the current experimental situation 
on the observed energy dependence for all the hadronic cross-sections including the 
photon induced ones. Then we will describe the original minijet model [7] which 
tries to calculate this dependence in a QCD based picture. After pointing out the 
problem of 'too fast' an energy rise predicted in these models, we will then show how 
the resummation of soft gluons can tame the rise [S] and how it is possible to obtain 
a satisfactory description of the current data with stable predictions at the LHC [S] 
within the framework. Then we turn to how the model can be extended to describe 
photon-induced processes in terms of the measured photon structure function and 
present our predictions for the high energy photon-proton cross-sections obtained in 
this model [6]. 

2 The data for pp to 77 

With the pp/pp , 7p and 77 cross-sections in the milibarn, microbarn and nanobarn 
range [TO], it is possible to accommodate them all in the same figure by multiplying 
the 7p and 77 cross-sections by 330 and (330) 2 respectively. This is an ad hoc factor, 
approximately given by Vector Meson Dominance (VMD) and the quark parton model 
[6] . Such compilation of all the proton and photon data for the hadronic cross-sections 
is shown in Figure HJ from [6]. While 7^ and pp/pp data could be interpreted as 
indicating the same rise, at least at presently reached accelerator energies, the figure 
seems to indicate that the a 77 data from LEP rise somewhat faster than the others. 
In the left panel of Fig. [T]the band corresponds to the predictions for the total pp/pp 
cross-sections according to the Block-Nordsieck (BN) model which we shall describe 
in the coming section. The same band is also shown in the right panel, together 
with predictions of other models [UJ for purely proton processes. The figure makes 
clear that any model has to address three issues: 1) what makes the cross-sections 
rise, 2) what makes them obey the Froissart bound and 3) whether they indicate any 
breakdown of factorisation between the proton and photon processes. 
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Figure 1: Proton and photon normalised cross-sections [6 J with the band of predic- 
tions for proton-proton from [9], also shown on the right, together with estimates 
from various models. 



3 Minijet models 

Minijet models were one of the early QCD based models which made an effort to 
understand the energy rise in terms of the rising gluon content of the hadrons at 
small values of x and the basic QCD cross-section [12]. The cross-section for the jet 
production in collisions of two hadrons A, B in the process: 

A + B ->■ X + jet (1) 

is obtained by convoluting the parton-parton subprocess cross-section with the given 
parton densities and integrating over all values of incoming parton momenta and 
outgoing parton transverse momentum p t , according to the expression 

a^(s,p tmin ) = dp t dx 1 dx 2 

J Ptmin 

x E MxuP^fM^p't)^ 1 (2) 

i,j,k,l 

where A and B are the colliding hadrons or photons. By construction, this cross- 
section depends on the particular parametrization of the parton densities evaluated 
at scale pf. This cross-section strongly depends on the lowest p t value on which one 
integrates, viz. p tm in- The term mini-jet was introduced long ago [33] to indicate all 
those low p t processes which are amenable to a perturbative QCD calculation but are 
actually not observed as hard jets; p t indicating the scale at which a s is evaluated in 
the mini-jet cross-section calculation. One can have ptmin ~ 1 -i- 2 GeV. This minijet 
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cross-section rises very fast with y/s, the rate of rise is controlled by the dependence 
of the parton densities at low-x values and ptmin used. In the simplest version of the 
model [7] a tot was assumed to be given by 

a tot = a + cT jet (s,p tmin ). 

Eventhough these early calculations caught the essence of QCD dynamics that can 
cause the rise of total and inelastic cross-section with energy, they predicted a rise with 
energy which depended on a parameter ptmin arbitrarily varying with energy. This 
was made necessary by the known low-x behaviour of the parton densities, which 
predicted a very fast rise with energy, leading to unitarity violation. The unitarity 
is achieved by embedding the minijet cross-section in an eikonal picture [14]. In fact 
many features of not just the cross-section but also quantities such as multiplicity 
distributions etc. can be successfully described in the eikonalised minijet picture [T5] . 

In the eikonal formulation, <j tot is obtained from an eikonal function x(b, s) which 
describes the impact parameter distribution during the collision, namely 

^,^/^b|l-e-(M |2 (3) 
atotai = 2 J d 2 b[l - e- Xmx ^cos^e X (b, s)} (4) 

CTtotal vaelasUc = duel = J <Pb[l ~ e" 2 ^^] (5) 

Neglecting the real part of the eikonal for the hadronic processes, one gets a very 
simple expression for the total cross-section 

dotal = 2 / d 2 b[l - (6) 

with n(b, s) the average number of inelastic collisions. 

Minijet models with the well motivated QCD input, embedded into the eikonal 
representation, are able to describe the early rise correctly. However, without any 
further energy dependent input they often fail to obtain the presently exhibited level- 
ling off at high energy, a behaviour already consistent with the Froissart bound [16] . 
It is clear that one needs additional input: within QCD there is one important effect 
which can modify the energy dependence, and this is soft gluon emission from the 
scattering partons. We shall discuss this effect in the next section. 

4 Block-Nordsiek Model 

The BN model differs from and improves on the usual mini-jet models, including the 
eikonalized ones, in three significant ways, namely 1) by implementing perturbative 
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QCD input through currently used PDFs for the mini-jet cross-sections, 2) intro- 
ducing soft gluon fc t -resummation to control the rise as the basic mechanism which 
describes the impact parameter distribution of the collision, with the upper scale in 
soft gluon resummation linked by kinematics to the mini-jet cross-section and finally 
3) pushing the soft gluon integral into the InfraRed (IR) region. We have discussed 
the mini-jet calculation in the previous section. Here we shall illustrate the basic 
features of resummation and our approach to it. 



4.1 Resummation 

We shall start by recalling some properties of soft photon resummation. In QED, the 
general expression for soft photon resummation in the energy-momentum variable 
can be obtained order by order in perturbation theory [T71 [TEJ HE] as 

d A P(K) = d 4 K J 0- e iK-x-h( X ,E) (7) 

where d A P(K) is the probability for an overall 4- momentum escaping detection, 

h(x, E)= f d 3 n{k)[l - e~ ik - x ] (8) 



and 

d 3 k 

d 3 n(k) = Wo \Uk,&})\ 2 (9) 

The electromagnetic current j M depends on the momenta of the emitting matter fields 
{pi} and to leading order in c^qed is given by 

«MMH-^I>^ do) 

with the sum running on all the entering particles (e, = +1) and antiparticles(-l). 
For outgoing fields, the signs are reversed. 

The expression in Eq.([7J can also be obtained employing the methods of statistical 
mechanics |20j . One starts with a discrete variable representation 

d'P{K) = Y,P{{n k })5\K -Y,n k k)d 4 K (11) 

■rik k 

and sums on all the possible Poisson distributed soft photon configurations, i.e. 

— n k 

P({n k }) = U k ^-exp[-n k } (12) 

where n k is the number of photons emitted with momentum k. From Eqs. f fTT|T2l 
one can determine the probability of observing a 4-momentum loss K accompanying 



5 



a charged particle reaction by using the integral representation for the 4-dimensional 
^-function and exchanging the order between forming the product in P({rik}) and the 
summation over all the distributions. Then, going from the sum over discrete values 
to an integral over soft photon momenta, one can obtain Eqs. (0 The derivation 
being semi-classical, it cannot give any information on d 3 n(k). For this, one needs 
to use the perturbative expression for the electromagnetic current as in Eqs. (|9][T0]) . 
Notice that in this derivation, it is energy momentum conservation which ensures the 
cancellation of the infrared divergence between soft and real quanta emission. This 
follows from the semi-classical approach, since in the infrared region the uncertainty 
principle does not allow to distinguish between real and virtual quanta. 
From Eq. (|TT|) . the integration over the 3-momentum variable gives 

dp M = ^TT«-<I> 5 (13) 

with 7 the Euler's constant, = (3({pi},m e ) defined as the integration over the soft 
photon angular distribution, i.e. such that 

/ d 3 n(k)=/3^- (14) 
Jn k k 

and E is the maximum energy allowed to soft photon emission. With E being the 
energy scale of the reaction, the integration of Eq. fll3j) up to the energy resolution 
AE, gives the well know behaviour (AE / E)- alogE ^ me proposed in the early days of 
QED [21]. 

While it is easy to obtain a closed form for the energy distribution, a closed form 
for the momentum distribution is not available. It is also not necessary, since the first 
order expression in oiqed is adequate. More interesting, for applications to strong 
interactions, is the transverse momentum distribution, namely 

d 2 P(K ± ) = d 2 K ± -±- [ d 2 h e -«-L-b-fc(MJ) (15 ) 
(27rr J 



with 



h(b, E)= J d 3 n(k) [1 - e iki b ] (16) 

For large transverse momentum values, by neglecting the second term and using a 
cut-off term as lower limit of integration, the above expression coincides with the 
Sudakov form factor [22]. In QCD, Eq. (JIB]) is used to discuss low-p* transverse 
momentum distributions through the expression 



where Cf = 4/3. 



, /, »^ 4C F f M n2 ,dk u 2M 

*<».*> =— U'^W^ (17) 
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We use resummation to approach the very large impact parameter region, which 
plays a role in such quantities as the total cross-section. For this, we need to explore 
the IR region and use the full range of integration of Eq. ( Tl6l) . In the region k t < A 
we propose to use an expression for the soft gluon spectrum which takes into account 
the effect of a rising confining potential. The expression to use will then be singular 
as kt —> 0, but must be integrable. As discussed in many of our publications [8], 
we shall interpolate between the asymptotic freedom region (AF) and the IR through 
the following expression for the strong coupling constant 

,, . 127T p , . 

^s(kf) = P 2 (18) 

( 33 -2A/)ln[l + p(^)p] 

which coincides with the usual one- loop formula for values of k t » A, while going 
to a singular limit for small kt, and generalizes Richardson's ansatz for the one gluon 
exchange potential to values of p < 1 [25J. However, notice that, in order to make 
the integral finite, we need p < 1. Also analyticity requires p > 1/2 [8J. 

Having thus extended the region of interest in Eq. (TTTf) to the large impact pa- 
rameter values, we use the Fourier transform of Eq. ( Fl5|) to describe the impact 
parameter distribution, and input it into the average number of hadronic collisions, 
n(b,s), generated by mini-jets, i.e. we write 

n har d{b, s) = A BN (b, s)a jet (19) 

with 

A BN (b, s) = A BN (b, M) = A e~ h ^ (20) 
with the normalization constant A 

A ° = 27rjbdb e~ h ^ M ) (21) 

and with the integral in Eq. ( 117]) extended down to k t = 0, to be used in Eq. (jSJ). 
The subscript BN indicates that this 6-distribution is obtained from soft gluon fc t - 
resummation into the IR region. 

Another important quantity in soft gluon resummation, is the upper limit of in- 
tegration, the scale E in QED calculation, and which we have indicated with M in 
the QCD integral. M is in general a function of the incoming and outgoing parton 
momenta, and can be determined by the kinematics, as discussed in [6] and following 
[26] . Upon integration over all the PDFs, it is seen to be of the order of ptmm and 
slowly varying with the cm. energy. 

4.2 BN model for protons 

In order to apply our QCD description, inclusive of mini-jets and soft fct-resummation, 
to scattering processes, we split the average number of collision n(b, s) as 

n(b, s) = n NP (b, s) + n hard (b, s) (22) 



7 



and parametrize the first term, which is meant to include all those processes which 
cannot be described through parton-parton scattering discussed in the previous sec- 
tion. At low energies, i.e. y/s approximately up to 20 -j- 30 GeV, n(b, s) n^fp(b, s), 
whereas as the cm. energy increases, the second term will asymptotically dominate. 
To see how, in the BN model, soft gluon resummation helps to tame the rise of the 
mini-jet cross-section and bring satisfaction of the Froissart bound, we shall thus look 
at the expression for a to tai at extremely high energies, namely 

a T (s) = 2 J d 2 b[l - e - nhard/2 ] (23) 

In the previous edition of this Conference, we have presented a preliminary version 
of the argument for satisfaction of the Froissart bound in our model [27]. A more 
accurate argument has now been given in [8]. Basically, we find that extending the 
integration in Eq. ( ITT}) into the IR region with our IR singular, but integrable, a s 
introduces a cut-off in 6-space, which behaves at least like an exponential (1 < 2p < 2). 
Such cut-off allows to evaluate the integral of Eq. ( 1231) at the value where the integrand 
suddenly decreases. Inserting the asymptotic expression for aj et at high energies, 
which grows like a power of s, and A BN (b, s) from Eq. ( 120]) . in such large-6, large- s 
limit, we obtain 



n hard = 2C(s)e-^ P (24) 



2/. 



where 2C(s) = A (s)ai(s/so) £ . The resulting expression for ot 



a T (s) « 2vr r db 2 [l - e - c ^- {bA)2P ] (25) 



o 



leads to 



with 



AV r (s) w (— ) / duu— = 2tiu /p (26) 
p Jo 



u = ln[^-}^elns (27) 
To leading terms in In s, we therefore derive the asymptotic energy dependence 

a T -)■ [eln(s)fM (28) 

Remembering that 1/2 < p < 1 [TJ, the above result shows that, with soft gluon 
momenta integrated into the IR region, k t < A, and a singular but integrable coupling 
to the quark current, our model leads to a behaviour consistent with the Froissart- 
Martin bound [lj. 
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4.3 The BN model for photons 



To apply the BN model to photon processes, we follow refs. (28J |29], and estimate the 
total cross-section as 

aZ = 2P had [ d 2 b[l - e- nlP ^' 2 ] , n»(b, s) = < P + A% N (b, s)%± (29) 

with Phad = 1/240 ~ 0{a) to represent the probability that a photon behaves like 
a hadron and with photon PDF's used to calculate the mini-jet cross-sections and 
the scale energy parameter M. A full discussion of how to evaluate this parameter , 
called q max when it is averaged over the densities, can be found in [6] . q max depends on 
the energy of the subprocesses and, being evaluated using the PDFs of the processes 
under consideration, depends on the specific choice of the parametrisation used for 
the parton densities in the photon and the proton. As q max increases with energy, the 
growth of the total cross-section due to mini jets is tempered by soft gluon emission. 
The calculated values of q ma x, for all the available parton densities reach some sort of 
saturation at high energies, which in turn reflects in the total cross-sections reaching 
a stable slope. 

The application of our model to photon total cross-sections shows some interesting 
features. While in the present accelerator energy range, ■yp data could be described 
also through factorization models, at very high energies in the TeV range, predictions 
differ. We show this in Fig. fl2]). The left panel plots the data up to the highest accel- 




Figure 2: yp total cross section from [6] in different energy ranges. 



erator energies. Data in this region come from cosmic rays [20] , from extrapolation of 
virtual photon data taken with the BPC (31] and from HI and Zeus experiments [32] . 
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In this energy range data can be accomodated by many models, including factoriza- 
tion models, in which models for proton cross-sections are extended to photons, by 
just multiplying the proton curve by a constant factor [33J or by extrapolating eikonal 
models with scaling factors in the impact parameter description of photons [31] or by 
assuming for photons the same rise with energy of proton cross-section [35] . However, 
as one can see from the right panel, the situation changes as the cm. energy of the jp 
system reaches into the TeV range. Here the discrepancy with factorization models, 
exemplified by the lower band, is clearly indicated both by the upper band obtained 
through our BN model, as well as by the dashed curve within it. This curve was 
independently obtained in [16] from a fit to accelerator data and confirms the validity 
of our model into an energy range so far inaccessible through particle accelerators. 

5 Conclusions 

We have discussed the results from a mini-jet model which incorporates soft gluon 
A; r resummation as a taming effect on the rapid rise with energy of low-x initiated 
mini-jet cross-sections. This has been applied to both pp/pp and 7p processes. We 
find that soft fc r resummation, inclusive of IR gluons with k t < A, plays a crucial 
role in trasforming the power like rise of the jet cross-sections into a more subdued 
logarithmic behaviour. We accomplish this through the use of a phenomenological 
ansatz for the coupling between soft gluons and the quark current which gives an 
expression singular but integrable. 
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